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Abstract 



Shell corrections of the finite deformed Woods-Saxon potential are calculated 
using the Green's function method and the generalized Strutinsky smoothing 
procedure. They are compared with the results of the standard prescription 
which are affected by the spurious contribution from the unphysical particle 
gas. In the new method, the shell correction approaches the exact limit pro- 
vided that the dimension of the single-particle (harmonic oscillator) basis is 
sufficiently large. For spherical potentials, the present method is faster than 
the exact one in which the contribution from the particle continuum states is 
explicitly calculated. For deformed potentials, the Green's function method 
offers a practical and reliable way of calculating shell corrections for weakly 
bound nuclei. 

PACS numbers: 21.10.Dr, 21.10.Ma, 21.60.-n 
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I. INTRODUCTION 



The positive-energy spectrum of the average single-particle potential plays a role in the 
description of weakly bound nuclei for which the Fermi level approaches zero (see Ref. 
For these nuclei, important for both nuclear structure and nuclear astrophysics studies, spe- 
cial care should be taken when dealing with the particle continuum which seriously impacts 
many nuclear properties, including bulk nuclear properties (e.g., masses, radii, shapes) as 
well as nuclear dynamics (i.e., excitation modes). 

In two earlier papers 0,0, a macroscopic-microscopic method was applied to nuclei far 
from the beta stability line. It has been demonstrated that the positive-energy single- 
particle spectrum does impact the results significantly, and the systematic error in binding 
energies, due to the neglect or the improper treatment of the particle continuum, can be 
as large as several MeV at the neutron drip line. In the first paper 0, both spherical and 
deformed nuclei were considered, and the continuum was approximated by a limited number 
of quasistationary states which resulted from a diagonalization of the Woods-Saxon average 
potential in a harmonic oscillator basis. In the vicinity of the neutron drip line, the result of 
the Strutinsky smoothing (standard averaging method) becomes unreliable and it deviates 
from the result of the semiclassical Wigner-Kirkwood expansion. The semiclassical method, 
which does not use the positive-energy spectrum explicitly, gives a more reliable estimate of 
the shell correction than the standard method (see Refs. [^Q and references quoted therein). 

In the following study ||^, carried out for spherically symmetric nuclei, a more detailed 
comparison was carried out between the Strutinsky smoothing method and the Wigner- 
Kirkwood expansion. In the Strutinsky method, the continuum effect was taken into account 
exactly, i.e., by calculating the continuum part of the level density from the derivative of 
the scattering phase shift with respect to single-particle energy. The smooth part of the 
continuum level density has been calculated by means of the contour integration along a 
path in the complex energy plane P]. Although the continuum level density was treated 
properly, it has been concluded that in most nuclei the plateau condition of the Strutinsky 
method could not be met. Therefore, this condition was replaced with the requirement 
of the linear energy dependence of the mean level density. This modification (which widens 
the range of the applicability of the Strutinsky procedure considerably) is referred to as the 
generalized Strutinsky method. The new procedure has proved to be very useful; in most 
cases it gives results reasonably close to the estimate of the semiclassical method. The 
exceptions are the neutron drip-line nuclei in which the neutron Fermi level approaches zero 
and the semiclassical procedure diverges 0. 

In the present work, the effect of the particle continuum on shell correction is further 
studied for both spherical and deformed nuclei by using the recently introduced Green's 
function approach The advantage of this method is that, by diagonalizing the finite 

single-particle potential in a square-integrable basis, one can get rid of the spurious contri- 
bution of the particle gas to the level density. We shall refer to this novel method as the new 
method in order to distinguish it from the commonly used standard smoothing procedure 
(or: old method) in which the spurious contribution from the particle gas is not subtracted 
(but diminished by using a reduced number of basis states). 

The paper is organized as follows. Section || contains a brief review of the shell-correction 
method and describes several versions of the smoothing procedure used. The numerical 
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results of shell-correction calculations are presented in Sec. |T| for both neutron-rich and 
proton-rich nuclei. Finally, conclusions are drawn in Sec. IV. 



II. STRUTINSKY SMOOTHING PROCEDURE 
A. Basic Definitions 



In the macroscopic-microscopic approach P,p|-p!^, the shell correction 



SEghell — Es,p, — -Es.p., (1) 

is defined as the difference between the total single-particle energy -Eg. p., 

^s.p. = E (2) 

i— occ 

and the smooth single-particle energy Es,p,. The shell correction represents the fluctuating 
part of the binding energy resulting from the quantal single-particle shell structure. 

For the sake of simplicity, we shall assume that the single-nucleon energy spectrum is 
that of a one-body Hamiltonian, 

H = f + V, (3) 

with an average single-particle potential V. In practice, the potential contains a deformed 
Woods-Saxon potential, the spin-orbit term and the Coulomb potential. Since the central 
potential is finite, the spectrum of H is composed of bound states with discrete negative 
eigenvalues (e, < 0) and the continuum of scattering states with positive energies (e > 0). 
Consequently, the single-particle level density is 

g{e) = fM{e)+g,{e), (4) 

where 

g^{e) = Y.2 6{e-e,) (5) 

i 

is the level density of the discrete (bound) states and g^e) is the continuum level density 
(it will be specified later). The factor 2 in Eq. (^ appears due to the two-fold Kramers 
degeneracy of the deformed single-particle energy levels. 

In the shell-correction method -Es.p. is calculated by employing the smoothed level 

density ^(e) obtained from g{e) by folding it with a smoothing function f{x): 

de' g{e')f( — ]. (6) 



7 J-oo \ 7 

In practical applications, for the folding function f{x) one is usually taking a product of a 
Gaussian weighting function, ;^exp(— x^), and a corresponding curvature correction poly- 
nomial of the order p |10| which is an associated Laguerre polynomial L^Jj^^{x) (p-even). The 
smoothed level density (Bl) defines both the smooth single-particle energy 
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Es.p. = / e~g{e)de, (7) 
and the smoothed Fermi level A. The latter is obtained from the particle number equation: 

N= ~g{e)de. (8) 



The smooth single-particle energy and the resulting shell correction have to be defined 
unambigously. Therefore, they must neither depend on the smoothing range 7 nor on the 
order p of the curvature correction. This requirement, referred to as the plateau condition, 
can be written as 

dE dE 

^^ = ,^^ = 0. (9) 

dj dp 

Of course, since one wants to eliminate the oscillations due to the shell structure, the smooth- 
ing range 7 should be greater than the average energy distance between neighboring major 
shells, huo^Al/A^/^ MeV 



For infinite potentials, such as an infinite square well, harmonic oscillator, and a deformed 
Nilsson potential, one can always find a range of the smoothing parameters 7 and p in which 
the smooth single-particle energy is independent of the values 7 and p |]14|,Q . For finite-depth 



potentials, however, additional complications arise due to (i) the presence of positive-energy 
continuum and (ii) the difficulties with meeting the plateau condition. We shall discuss 
these points in the following. 



B. Effect of the Unbound Spectrum 

The need for calculating the continuum level density, gd^), appears whenever one deals 
with finite-depth potentials. For spherically symmetric potentials, the continuum level den- 
sity H,p!5Hl8| is defined by means of the scattering phase shifts Sij{e): 

=ii:(2j+i)^- (10) 

vr , . ae 

'J 

For realistic nuclear potentials, phase shifts have to be calculated by numerically solving the 
radial Schroedinger equations and by matching the wave function to the asymptotic solution 
at a distance where the nuclear potential can be neglected. This procedure has to be carried 
out for every partial wave below a certain angular momentum cut-off on a fine mesh in the 
positive-energy region. In order to prevent sudden jumps in gd^) around narrow resonances, 
a new calculational method employing the Cauchy theorem was introduced in Ref. 0. Here, 
the complex energies Wi of the Gamow resonances (poles of the S-matrix) are localized first, 
then a contour of the complex energy plane is chosen. The contour, denoted by L, should go 
far away from the poles. The mean level density g{e) is then calculated as a sum over bound 
and those resonant states which lie between L and the real energy axis and an integral term 
along a contour: 
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m = E / j + X 3c{w)f l^^j . (11) 

Apart from the numerical errors, this procedure gives the continuum level density exactly. 
Therefore, we call this approach as numerically exact or, simply, exact. 

For deformed single-particle potentials, the continuum level density has a more compli- 
cated form and can be expressed by the on-shell S-matrix S{e, k, k') as 

9c{e) = ^Tr[^(e, k, fc')*^5(e, k, k')]. (12) 

If one wants to use expression (^) for calculating the continuum level density, one has to 
determine the S-matrix by solving the coupled system of differential equations for each value 
of e 10] . In practice this is a difficult task. 

In the old method, the single-particle Hamiltonian is diagonalized in a square-integrable 
basis formed from the eigenstates of an infinite potential. This potential can be either a finite- 
range potential contained in an impenetrable box (i.e., an infinite wall at a certain distance) 
or the harmonic oscillator potential. Since the number of basis states is always assumed 
to be finite, the diagonalization results in a discrete set of eigenstates. The eigenstates 
with negative energy approximate the bound states of the original Hamiltonian, while the 
positive-energy quasi-bound states mock-up the effect of the particle continuum in a very 
crude way. It has early been realized that in the application of this method one should not 
use too large a basis; otherwise, the level density around the zero energy would increase 
dramatically. (In fact it diverges as the basis size goes to infinity.) In order to avoid this 
catastrophe, the use of a harmonic oscillator basis with 12-14 harmonic oscillator shells was 



recommended [|rT|]. One of the objectives of this paper is to perform the critical evaluation 
of the standard smoothing method by calculating the continuum level density in a more 
reliable way using the Green's function approach. 



C. Green's Function Method 



Although the Green's function approach to the single-particle level density was developed 
long ago (see, e.g., Refs. |[T5|JT9| -^ ) , it is somehow surprising that so far it has not been 
widely applied. Below, we briefiy summarize the main features of this method. More details 
can be found in Refs. 0J^. 

A Hamilton operator with an infinite potential, Hoo, has only discrete energy eigenvalues 
and its eigenf unctions are all square integrable. Therefore, in this case, the single-particle 
level density is given by Eq. (^). By introducing the Green's operator, Goo{z) = (2 — i^oo)~^ 
(7d(e) can be written as 



9d[e) 



-ImjTr [G'oo(e)l}. 



(13) 



As discussed in Ref. 
density (^ becomes 



for a Hamiltonian H containing a finite potential the full level 



— Im |Tr 



G+(e) - 



(14) 
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where G~^{z) = {z — H + iO)~^ and Gq{z) = {z — Hq + iO)~^ is the free outgoing Green's 
operator associated with Hq = T. The interpretation of Eq. ( |I^ ) is straightforward: the 
second term contains the contribution to the single-particle level density originating from 
the gas of free particles. 

Let us now introduce an approximation to the exact expression ([M])- To this end, we 
diagonalize H and Hq in an orthonormal basis formed from the M square-integrable basis 
functions. The resulting approximate eigenenergies of H and Hq are denoted by and e°, 
respectively (z=l,...,M). This procedure amounts to a projection of both Hamiltonians into 
the M-dimensional Hilbert space of square-integrable basis functions. The level density 
can then be approximated by the difference of the discrete level densities of the two projected 
Hamiltonians: 

M M 

guie) = Y: 25(e - e.) - J] 25(e - e°). (15) 

i=l i=l 

By increasing the dimension M, the bound eigenvalues of H converge to the exact single- 
particle energies while the positive-energy eigenvalues will tend to approach zero energy. 
The eigenvalues e^, which are obviously different from the positive energies e, for any finite 
M, can, in fact, compensate for the spurious increase of the level density around the zero 
energy if the smoothing procedure (H) is applied: 

~9m{^) = - de' gMie')f ('^] . (16) 

It has been shown in Ref. |^ that the exact smoothed level density g{e) can be reproduced 
by gui^) in the limit of large M. 

In this work we employ the (stretched) harmonic oscillator basis. As discussed below, 
thirty oscillator shells are sufficient to guarantee the convergence of results. In the standard 
(old) method, one takes much fewer states (12-14 oscillator shells) and the second term in 
Eq.(p!5D is not subtracted. Consequently, the results are not stable as M is varied. 



D. Generalized Shell-Correction Method 

In the standard Strutinsky smoothing method, when applied to finite potentials, it is 
difficult to meet the plateau condition (^ [g]. The more detailed study of Ref. demon- 
strated that the plateau condition can seldom be satisfied even if the particle continuum is 
properly accounted for. Fortunately, it is possible to replace the standard plateau condition 
with a new requirement which yields unambigous shell-correction values [Q. By comparing 
the smoothed Strutinsky level densities with those obtained in the semiclassical Wigner- 
Kirkwood method, it was found that they are in good agreement, apart from the low and 
the high ends of the spectra. In the intermediate energy region, the average level density 
shows linear dependence on e. (The linearity of the semiclassical level density for heavy 
nuclei was noticed already in Ref. |jl5|.) Guided by this observation, the shell-correction 
method was generalized by replacing the plateau condition with the requirement that in an 
energy interval [e/, £„] which is wider than the average distance between neighboring major 
shells. 
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e„ - e; = 1.5 %uq, (17) 

the deviation of g{e,'y,p) from hnearity should be minimal [Q. In practice, one has to 
minimize the deviation 

X^{l,p)= [g{^,l,P) - a-bef de, (18) 

where the parameters a and b are uniquely determined for each value of 7 and p by the 
method of least squares. 

Figure |l] displays ^ function of 7 and p for the neutrons in the spherical superheavy 

nucleus Z=114, iV=184. It is seen that has two minima for each p value, and there is a 
clear correlation between p and 7 (minima are shifted to larger values of 7 with increasing 
p). The position of the first minimum is between 1.12 and 1.48 Hujq (i.e., between 6.6 MeV 
and 8.6 MeV). The second minimum appears at larger 7 values, namely between 7.7 MeV 
and 12 MeV. As demonstrated in Ref. 0, the average level densities corresponding to the 
first minimum of practically do not depend on p in the negative energy region. This is 
also valid for the second minimum in x^. However, a difference can be seen if one compares 
average level densities calculated at different minima. The inset of Fig. |I] shows ^(e) for 
both minima at p=10. The level density corresponding to the lower-7 minimum preserves 
its linearity for a wider range of energies, and the linearity is best fulfilled in the energy 
region which is midway at the bottom and the top of the potential. Therefore, in our 
calculations [Eq. ([TsD j we fixed the energy interval so that it is centered around the half of 
the energy of the lowest single-particle level. Far from this central region g{e) varies rapidly 
and this forces the smoothed level density to oscillate. These oscillations can be viewed as 
end effects, and they are largely independent of the shell structure. For example, for the 
harmonic oscillator potential whose spectrum has no natural upper bound, these oscillations 
occur around the bottom of the potential well. For the Woods-Saxon potential, additional 
oscillations occur around the threshold energy. 

It is well known that the realistic value of the smoothing parameter 7 has to lie in a 
certain energy interval The value of 7 should be large enough to wipe out shell effects 



in the energy range of a typical distance between shells, but it should not be much larger 
to avoid bringing the threshold oscillations down to lower energies. To prevent this, in our 
calculations we always use the 7 values corresponding to the first minimum of x^. For the 
case shown in Fig. |I|, shell correction changes by 0.3 MeV if one uses the value of 7 at the 
second minimum instead. While in this case the change in shell correction is well within the 
uncertainty of the model, the difference is more pronounced for lighter nuclei. For the very 
light nuclei, the end effects dominate the energy dependence of ^(e) in the whole energy 
range. Consequently, the Strutinsky smoothing method cannot be meaningfully applied to 
these systems. 

III. RESULTS 
A. Model Parameters 

In the calculations presented in this work, we have used the average, axially deformed 



Woods-Saxon (WS) potential of Ref. ||2^, which contains the central part, the spin-orbit 
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term, and the Coulomb potential for protons. The potential depends on a set of deformation 
parameters, /3a, defining the nuclear surface: 



R{e;(3) = C{f3)roA'/' 



l + E/5Ano(^) 

A 



(19) 



where the coefficient C assures that the total volume enclosed by the surface (|T9|) is con- 
served. The Coulomb potential has been assumed to be that of the charge {Z — l)e uniformly 
distributed within the deformed nuclear surface. We employed the set of WS parameters 
introduced in Ref. For details pertaining to the WS model, see Refs. |p|,^. 



The deformed WS Hamiltonian was diagonalized in the deformed harmonic oscillator 



basis using the computer code of Ref. ||23|. For the diagonalization we took all the (stretched) 
oscillator states having the principal quantum number less or equal than A/'max- (In short, 
we took A^max deformed shells.) The diagonalization of Hq was carried out in precisely the 
same basis. 

When adopting the new scheme, the important question is how many harmonic oscillator 
shells are needed in order to reproduce the exact results of Ref. 0. Naturally, the value of 
Nraax depcuds ou the size and the shape of the potential to be diagonalized, and also on the 
oscillator frequency hu = rjTiujQ (the default value of rj is 1.2). 

The convergence of 5i?sheii for neutrons as a function of N^^^k is illustrated in Fig. |^ for 
^^^Sn and ^^^Sn (which is weakly bound). One can see that for both nuclei the shell correction 
obtained in the new procedure quickly converges to the exact value, and at A^max=30 the 
agreement is very satisfactory. Therefore, in the following, we shall use 30 oscillator shells 
when applying the new method. 

As expected, the results of calculations using the standard method do not stabilize 
with A/'max- However, for the recommended values of A^max=12 and 14, the shell correction 
produced with the old method differs from the exact value by less than 1 MeV. However 
this apparent agreement seems to be accidental. We display in Fig. |^ and ^d, respectively, 
the single-particle energy and the smoothed single-particle energy for ^^^Sn as a function 
of A^max- One can notice that at A^max=12 the total single-particle energy differs from the 
exact value by about 8 MeV. Since the corresponding smoothed single-particle energy is 
also shifted by about 8.5 MeV, the resulting shell correction differs only by 0.8 MeV from 
the exact value. Consequently, an acceptable agreement for 5-Esheii comes as a result of 
cancellation between two large numbers, each subject to large errors. As one approaches 
the neutron drip line, the accurate calculation of single-particle energies requires a rather 
high number of shells and/or a basis optimization with respect to the parameter rj which 
determines the oscillator length. This is illustrated in Fig. § which shows the convergence of 
the total neutron single-particle energy for ^^°Zr at a large quadrupole deformation /32=0.6. 
Clearly, for a weakly bound and deformed nucleus one needs at least A^max=30 oscillator 
shells to reach the convergence with the standard value of r]. Of course, by increasing the 
oscillator length, i.e., by choosing a smaller value of r], one can improve the convergence 
significantly for a system with a spatially extended density. In this example, one can arrive 
at a reasonably accurate value of E^.-p. by using A^max=20 and ?7=0.8. Another practical way 
of improving the convergence of single-particle energies is to use the modified oscillator basis 
obtained by means of the local scaling transformation [p5|] . 

Figure § shows the neutron smoothed level density for ^'^^Sn, a relatively well bound 
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nucleus, calculated with different methods. The new method with A^max=30 describes very 
well g{e) in the whole region of negative energies. This proves that the Green's function 
approach can be used with confidence, even for weakly bound systems. On the other hand, 
the average level density obtained with the standard method never stops increasing, and its 
deviation from the exact result shows up already at e=-18 MeV. The result displayed in 
Fig. ^ demonstrates that even for well-bound nuclei, the shell corrections calculated using 
the old method are prone to significant errors. 

B. Deformation Effects 

In order to investigate the deformation dependence of shell corrections, we performed 
calculations for ioO'iio,i202j^ function of (32 (other deformation parameters were assumed 
to be zero). In the Green's function variant, the generalized plateau condition was used; the 
resulting values of 7 were also employed in the standard Strutinsky calculations. The results 
are shown in Fig. ^ As expected, the most pronounced difference between the results of the 
two methods is for the weakly bound nucleus, ^^°Zr. This difference does depend on P2', part 
of it can be attributed to the deformation dependence of the smoothing width. (It should 
be noted that in the deformed calculations of Ref. [0 7 was assumed to be constant.) 

Since the effect of the particle continuum on SEgheii should be less pronounced for systems 
that are bound better, one would expect the two methods to yield similar results for lighter 
Zr isotopes. However, as seen in Fig. ^, the difference between both methods is negligible 
only at very low energies, e<18MeV. For higher values of e (or Fermi level A), the difference 
between the smoothed level densities is not negligible and it is not even a monotonous 
function of A. As a result, one can notice in Fig. ^ that for ^^°Zr the results of both 
methods are very close while in a better bound nucleus of ^''^Zr they differ more. It is 
interesting to note that for ^''^Zr and ^^°Zr the deformation dependence of SEsi^eii is very 
similar in both variants of calculations. This is in accord with the observation made in 
Ref. 10] that the difference between the shell corrections obtained in the standard method 
and the semiclassical approach depends rather weakly on deformation. 

As discussed earlier, two main sources of the error of the standard smoothing procedure 
are (i) the error in determining the total single-particle energy, and (ii) the uncertainty of 
smoothing that influences the value of the smooth single-particle energy. Since (i) and (ii) 
are not independent, a large cancellation takes place which might reduce the total error. 
However, it is a priori difficult to predict how large this difference is and how strongly it 
would affect the predicted position of the drip line. In order to shed some light on this 
question, we carried out a comparison between shell corrections calculated with the two 
methods for the spherical Sn and Z=114 isotopes, and for the deformed Zr and Er isotopes. 
(In the latter case, we fixed deformation at /32=0.2.) These nuclei represented medium- 
mass and heavy nuclei where the generalized shell-correction method can be applied. The 
results are presented in Fig. |^. Except for well-bound Zr and Z=114 isotopes, the difference 
between the shell corrections calculated using the old and the new methods are on the order 
of MeV, and, except for the Sn isotopes, it is rather large when approaching A = (drip line). 
Although our calculations do not aim at determining the actual position of the drip line, they 
give a reasonably good estimate for the uncertainty of the old procedure. If one identifies 
the drip line with the neutron number where A becomes positive (this assumption is usually 
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violated in actual calculations because of the lack of self-consistency between the microscopic 
and macroscopic parts of the energy formula 0]), our limited calculations suggest that the 
drip-line predictions by the standard method are prone to severe uncontrolled uncertainties. 

C. Modification of the Green's Function Method to the Proton Case 

In the presence of the long-range Coulomb potential, the free Hamiltonian appearing 
in Eq. ( p!4D has to be modified. Indeed, for the protons, the asymptotic behavior of the 
scattering states is that of the Coulomb functions, not plane waves. Therefore, in this case, 
for the free Hamiltonian we take 

Ho = f + Vcouu (20) 

with Vcoui being the Coulomb potential. The role of the Coulomb term is to effectively push 
the continuum up in energy to the top of the Coulomb barrier. The results are insensitive 
to the radius of the Coulomb potential in the free Hamiltonian. As a matter of fact, even a 
point Coulomb potential can be used in Eq. (pO|) [p6| . 

Figure |^ shows the smoothed proton level density in the proton-rich nucleus ^^''Pb. One 
can see that already with a rather low value of iVmax=19 the new method reproduces the 
exact smoothed level density in the whole range of negative energies, and A^max=30 gives an 
excellent agreement with the exact result. The reason that relatively low values of A^max are 
sufficient in the proton case is that even slightly unbound states (narrow proton resonances) 
are well localized due to the confining effect of the Coulomb barrier. 

IV. CONCLUSIONS 

In this work, we employed the Green's function oscillator expansion method to calcula- 
tions of shell corrections. For spherical nuclei, the new method has proved to be a fast and 
very accurate approximation to the exact procedure. It also allows for a straightforward 
generalization to deformed shapes. In essence, the method is based on two simultaneous 
diagonalizations in a large oscillator basis. The first diagonalization involves the actual one- 
body Hamiltonian while the other one is carried out for the free Hamiltonian representing 
the particle gas whose contribution to the level density should be subtracted. For the neu- 
trons, the free Hamiltonian is given by the kinetic energy operator, while for the protons it 
also includes the Coulomb potential. In practice, the space of 30 (stretched) oscillator shells 
is sufficient to guarantee the stability of results. This relatively large (but still tractable) 
space is necessary not only for the proper treatment of the free gas but also for the accurate 
calculations the total single-particle energy. 

As demonstrated in our study, the use of the standard smoothing procedure can lead 
to serious deviations when extrapolating off beta stability. In particular, the particle drip 
lines predicted in the traditional approach can be very uncertain. (The systematic error in 
6Eshe\\, due to the particle continuum, can be as large as several MeV at the neutron drip 
line.) According to our calculations, the error on dEg^eW depends weakly on deformation in 
most cases. It is only for the weakly bound nuclei that the difference between the old and 
new methods exhibits a sizable deformation dependence. 
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There is no simple "fix" that would cure the deficiencies of the standard Strutinsky 
procedure when applied to finite-depth potentials. One does need the large basis in order to 
guarantee the stability of -Eg. p.- On the other hand, at these large values of A^^ax, the smooth 
single-particle energy becomes unreliable due to the unphysical increase of the quasi-bound 
levels around the threshold. We believe that the new Green's function method, together 
with the generalized plateau condition, is a very useful tool that should be employed in 
future global calculations of nuclear masses in the framework of the one-body (macroscopic- 
microscopic) description and in level-density calculations for spherical and deformed nuclei. 
Of course, the new procedure does not remove the generic problem of the lack of the self- 
consistency condition between the microscopic and macroscopic Fermi energies 0. Recently, 
the Green's function method, based on self-consistent potentials obtained in Hartree-Fock 
and relativistic mean-field calculations, was used to extract shell corrections in the spherical 
superheavy nuclei . Although these calculations were not done using the oscillator basis 
expansion method but directly in the coordinate space, their main principle is the same as 
that discussed in this paper. 
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FIG. 1. Dependence of of Eq. ( |T^ ) on the smoothing range 7 for the spherical superheavy 
nucleus Z=114, A'^=184. The calculations were performed for several values of the curvature order 
p. The inset shows the average level densities corresponding to the two minima A and B of with 
p=W. 
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FIG. 2. Dependence of the neutron shell correction on the size of the harmonic oscillator basis 
for (a) ^'^^Sn and for (b) ^^^Sn. All of the single-particle states with principal oscillator quantum 
number less or equal than A'^max were considered in the diagonalization. The Amax-dependence of 
-Eg. p. and E^.p. for ^^^Sn is shown in portions (c) and (d), respectively. 
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FIG. 3. Dependence of the total single-particle energy for i^"Zr at p2=0.6 on the size of the 
stretched harmonic oscillator basis A^max for the two values of the oscillator frequency parameter 
{hu/nuJo=0.8 and 1.2, where Huq = 4l/A^/^ MeV). 
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FIG. 4. Neutron average level density for the spherical nucleus ^^^Sn calculated using the 
Green's function method with A'inax=30 (dashed line), the standard smoothing method with 
-/Vmax=12 (dotted line), and the Gamow-state technique (exact result, solid line). 
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FIG. 5. Neutron shell corrections for ^'^'^Zr (short dashed lines), ^^^Zr (dashed lines), and ^^''Zr 
(solid lines). The results of the Green's function method with iVmax=30 are shown by thick lines 
while those of the standard smoothing method with A^max=12 are indicated by thin lines. 
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FIG. 6. Neutron shell corrections for spherical Sn and Z=114 nuclei and deformed (/52=0.2) 
Zr and Er isotopes as a function of the neutron Fermi level A. The Green's function and standard 
calculations are shown by filled and open symbols, respectively. 
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FIG. 7. Proton average level density for the spherical nucleus Pb calculated using the Green's 
function method with A^max=30 (dashed line) and A'max=19 (dotted line), and the Gamow-state 
technique (exact result, solid line). 
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